Gated Silicene as a tunable source of nearly 100% spin-polarized 

electrons 



2 



Wei-Feng Tsai, 1 Cheng- Yi Huang, 1 Tay-Rong Chang, 
Hsin Lin*, 3 Horng-Tay Jeng, 2 ' 4 and A. Bansil 3 

1 Department of Physics, National Sun Yat-sen University, Kaohsiung 80424, Taiwan 
2 Department of Physics, National Tsing Hua University, Hsinchu 30013, Taiwan 
3 Department of Physics, Northeastern University, 
Boston, Massachusetts 02115, USA 
Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 
(Dated: March 13, 2013) 

Abstract 

Silicene is a one-atom-thick 2D crystal of silicon with a hexagonal lattice structure 
that is related to that of graphene but with atomic bonds that are buckled rather than 
flat. This buckling confers advantages on silicene over graphene, because it should, in 
principle, generate both a band gap and polarized spin-states that can be controlled 
with a perpendicular electric field. Here we use first-principles calculations to show 
that field-gated silicene possesses two gapped Dirac cones exhibiting nearly 100% spin- 
polarization, situated at the corners of the Brillouin zone. Using this fact, we propose 
a design for a silicene-based spin-filter that should enable the spin-polarization of 
an output current to be switched electrically, without switching external magnetic 
fields. Our quantum transport calculations indicate that the proposed designs will be 
highly efficient (nearly 100% spin polarization) and robust against weak disorder and 
edge imperfections. We also propose a Y-shaped spin/ valley separator that produces 
spin-polarized current at two output terminals with opposite spins. 
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I. INTRODUCTION 



Controllable non-magnetic spintronic devices are highly desirable in the development of 
spintronics applications which is an area of intense current interest [l- 3)]. While the spin 
polarization is generated by spin-orbit coupling (SOC), the entanglement between spin and 
orbital degrees of freedom due to the SOC significantly reduces the degree of spin polarization 
of spin-split states in most non-magnetic semiconductors including the topological insulators 
41. Recently, silicene, a close relative of graphene {s|, has been predicted and synthesized {§- 
16| . The band structure of silicene is similar to that of graphene in that the conduction and 
valence edges occur at the corners (K and K' points) of the Brillouin zone. Silicene however 
possesses fundamental differences fro m g raphene driven by the presence of a larger SOC, 



which opens gaps at the K-points 17|, [l8|. These gaps can be tuned with an external E-field 



perpendicular to the plane, which breaks the inversion symmetry (IS) of the system due to 



the presence of buckling in the honeycomb structure 19| . In this way, silicene can overcome 



difficulties associated with graphene in electronics applications (lack of a controllable gap) 
2ol 2lj |. potential applications of graphene in nano-electronics due to the available spin, 



valley and pseudo-spin degrees of freedom notwithstanding [22|-|27l|. 

In this work we present detailed first-principles computations to show that the band 
structure of gated silicene harbors two nearly 100% spin-polarized Dirac cones at the K- 
points. We use this key finding to propose designs of a silicene-based spin-filter as well 
as a spin-separator. Our quantum transport calculations show how a silicene based high 
efficiency spin filter (nearly 100% spin polarization) suitable for high-frequency electronics 
applications can be designed to switch the output spin current simply by gating without the 
need to switch magnetic domains. We analyze reflection and transmission of current at the 
junction of two domains with different external fields in order to obtain insight into trans- 
port selectivity between spin and valley degrees of freedom. High efficiency and tunability of 
our spin filter takes advantage of bulk charge carriers rather than the edge current in quan- 
tum spin Hall (QSH) systems, making our silicene filter robust against weak disorder and 
edge imperfections. We also propose a Y-shaped spin/valley separator that produces spin- 
polarized current at two output terminals with opposite spins. This three terminal device 
could be used for logical circuits beyond binary operations. Ge, Sn, and Pb counterparts 
of silicene are shown to have similar properties, but their larger SOC results in larger en- 



ergy differences between the spin-split states making these materials better suited for room 
temperature applications. Silicene, germanene, Sn and Pb thin films would thus be ideal 
materials for functional electronics and spintronics applications. Silicene however stands out 
because it could be integrated more naturally into the existing Si-based technologies. 



II. RESULTS 



Band structure and phase transition via IS breaking. Fig. [T] shows the low- 
buckled honeycomb lattice of silicene with its two sublattices A and B, and key features of 
the band structure. At low energies, the 7r-electrons dictate the system and reside around 
two inequivalent valleys, K and K' points, in the first Brillouin zone. Although we have 
carried out first-principles computations, insight into the underlying physics is obtained by 



describing the system by a minimal, four-band tight-binding model 171. Il8] 
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where c ia creates an electron at site i with spin polarization a. The first term is a nearest- 
neighbor (NN) hopping term on the honeycomb lattice. The second and third terms are 
intrinsic (IS preserved) and Rashba (IS preserved but z — > —z mirror symmetry violated) 
spin-orbit interactions, respectively, which involve spin dependent next nearest-neighbor 
(NNN) hopping, s are the Pauli spin-matrices, = (dj x dj) z /|dj x dj| = ±1 with two 
NN bonds, dj and dj, connecting the NNN d^. //^ = 1(— 1) when connecting A-A (B-B) 
sites. Notably, these SOCs, originating from buckling of the structure, are what distinguish 
silicene from graphene, even though both materials form a honeycomb lattice. The fourth 
term is an IS broken, staggered sublattice potential (£j = ±1), which arises when an external 
out-of-plane electric field, E z , is applied. The fifth term represents the effect of an applied 
out-of plane Zeeman (exchange) field with strength h. The last two terms, corresponding 
to the microscopic responses to the external fields, turn out to be the driving forces for the 
phase transition discussed below. A phase diagram is provided as Supplementary Fig. SI in 
the Supplementary Information. All the coupling parameters (Table I) have been obtained 



by fitting to first-principles band structures of stand-alone silicene sheet and Ge, Sn, Pb 
counterparts. Note that the fourth term in our Hamiltonian of Eq. ([!]) refers to the net 
electrostatic potential difference due to the environment, and accounts for effects of both 
the external field and the presence of the substrate. We expect therefore our modeling to 
be relevant for silicene placed on a substrate so long as the coupling with the substrate is 
not so strong as to modify the electronic structure of the silicene sheet. Although the band 
gap could be enhanced through strong interaction with the substrate, which is desirable 
for room temperature applications, it will be a tradeoff between achieving a larger gap and 
the possible reduction in the degree of spin-polarization resulting from the involvement of 
substrate orbitals. Moreover, the origin of the band gap itself could be different for different 
substrates. In order to help guide search for a suitable substrate, Supplementary Fig. S2 
gives the size of the band gap as a function of the lattice constant. It would perhaps be best 
to proceed by placing the silicene sheet on an insulator or semiconductor substrate where the 
lattice mismatch is not too large. It will be interesting however to further examine strong 
substrate coupling effects on a silicene sheet. Finally, we note that if a non- magnetic metal 
with small spin-orbit coupling is used as contact material, the spin-up as well as spin-down 
states of silicene will couple with more or less equal strength with the substrate and as 
a result spin-transport will not be affected much even though contact resistance may be 
altered substantially. 

In the absence of external fields, two of the four bands are occupied in Eq. (JTJ, so that 



the ground state of silicene is a QSH insulator with a SOC gap of 2Aso=8.4meV 



13, 



19]. 



The topological (Z2) nature of this ground state can be easily examined by either the parity 



analysis 30| or by showing the presence of helical edge states in a zigzag strip geometry. Note 
that the spin degrees of freedom in the band structure are still degenerate, as a consequence 
of both time reversal symmetry (TRS) and IS, even though is no longer a good quantum 
number in the presence of NNN Rashba SOC (except at K-points where Ar term makes no 
contribution) . 

When E z is applied to the buckled structure in which the two atoms within the unit 
cell are not coplanar, we naturally obtain a non-vanishing A v = IE Z caused by IS breaking. 
The estimated value of I is given in Table I. At K point, the energy gap, now 2|Aso — A v |, 
decreases linearly as E z increases until a critical field E c is reached. More importantly, the 
spin degeneracy is lifted with energy splitting of 2A V as shown schematically in Fig. [lb and 



precisely in Fig. le. We emphasize that this splitting arises from the spatial IS breaking as 
well as the presence of SOC in silicene. The spin polarization profile around the K' point is 
opposite to that of the K point as required by TRS. Furthermore, by symmetry arguments 
mentioned above, remains a good quantum number at K-points. 

As approaching the critical field, E z = E c , the gap shrinks to zero, forming a Dirac- 
like cone near each valley with the spectrum, y/k 2 + m 2 , where m = a/ A 2 , + a 2 \\k 2 and k 
denotes the momentum with respect to the K-point (See Fig. [lb). In sharp contrast to the 
Dirac cone in graphene, where each cone is spin degenerate, here the spin is fully polarized 
along the z direction at each valley (Fig. Id) and develops a small in-plane component of 
the form, a(s x k) 2 , with a of order Ar. Therefore, in the critical phase (dubbed spin-valley- 
polarization metal (SVPM)[3l|) near half- filling, the physics of the system is dominated by 
two nearly fully spin-polarized (with opposite polarizations) Dirac-like cones at K and K' 
points. 

For E z > E c , the energy gap re-opens at each valley and drives the QSH phase into 
the topologically trivial band insulating (BI) phase. It is interesting to note that the spin 
polarization is unchanged around K-point, although the probability amplitudes for A and B 
sublattices in the wavefunction is different. Starting from the BI phase with spin polarization 
of each band similar to that shown in Fig. [Tb, a Zeeman field h along -z direction cannot 
affect the direct band gap at each valley. However, the conduction band bottom would shift 
downward at K' point, while the valence band top would shift upward at K point. Increasing 
h therefore reduces the indirect band gap of the system and drives it into a valley-imbalanced 



metallic state, dubbed valley-polarized metal ( VPM) 3jJ , after crossing the zero indirect gap 
point. 

Field-tunable spin filter. Because various phases in silicene can be realized through 
the interplay of applied electric and magnetic fields 31], it becomes possible to set up a 
silicene based device working as a high efficiency spin filter. We illustrate this possibility by 
considering a 2D device consisting of a quantum point contact (QPC) in a silicene thin film. 
As seen clearly in Fig. [2^, a QPC is characterized by a short and narrow constriction. When 
a current / passes through the device by the application of an applied potential across the 
two wide regions, the conductance is sharply quantized in units of e 2 /h. To make the valley 
degrees of freedom well separated and thus keep the features unique, we choose the zigzag 
edges for the whole geometry along the direction of current flow. Note that a similar device 



has been proposed to realize a valley filter in graphene, but that works by using special 



23]. 



properties of the edge state and it is completely different from our work here 

We now show through quantum transport calculations that on the first few conductance 
plateaus the QPC produces an almost fully spin-polarized current. By locally changing the 
potential barrier via gating control in the constriction, the spin polarization direction can 
be easily reversed. To begin with, we map the tight-binding description of silicene, Eq. ([Q), 
to the geometry shown in Fig. EK The two opposite wide regions in the SVPM phase are 
arranged with Fermi energy Ep = 0.07t to model metallic source and drain (Ep = defines 
position of the "Dirac point"). Note that in the source region both spin polarizations are 
present (See Fig. [2b). In the constriction, we keep the system in the marginal VPM phase 
with a fixed E z > E c and the applied Zeeman field h = A v — Aso such that the indirect 
gap becomes zero (See Fig. [2hi) . Furthermore, we add an electrostatic potential barrier 
U(x{), which varies only along the current flow direction and is non- vanishing only in the 
constriction region, 

U( Xi ) = U [Q Ls {xi) - e Ls (x, - L)], (2) 



where L is the length of the narrow region, L s is a parameter for smoothing the potential, 
Q Ls ( Xi ) = (1) for ^ < -LJ2 (a* > L s /2) and Q Ls (xi) = § + §sin(^) for \ Xi \ < L s /2 ( Xi = 
\ defines the middle of the constriction). Two types of potential barriers are considered: 1) 
L s = 0, a rectangular shape with potential height Uq and 2) < L s < L, a smooth shape, 
as shown in Fig. 2a (lower part) by solid and dashed curves respectively. 

A typical dispersion relation for the wide region in the SVPM phase is explicitly shown 
in Fig. [2b. Given E-p, this region provides a total of 2N right-moving modes. In particular, 
as long as E F < (A v + Aso) an d Xr <C {Aso^K N propagating modes coming from the left 
valley will carry only down-spins and the other half modes coming from the right valley will 
carry only up-spins. This result has the merit of conveniently counting the spin polarization 
defined below, though, not necessary in our calculation. Note also that the two states closest 
to k = it I a have strong edge state character. On the other hand, for the constriction the 
dispersion relation is shown in Fig. EH. In the marginal VPM phase, each valley still has a 
direct gap but is shifted upward or downward due to the presence of h, resulting in a zero 
indirect gap. The two-terminal conductance of the QCP can be calculated by the Landauer 
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formula 

G = ^£ivi 2 ^(T t +n), (3) 

flU 

where the spin- resolved transmission probability T-^q = 5^metQ0 ^« l^mn| 2 with m and n 
representing outgoing and incoming channels, respectively. The transmission matrix t mn 
can be computed numerically by the iterative Green's function method (32}] . The spin 
polarization can now be expressed as 

P=^-^- (4) 

For < P < 1, the transmitted current is polarized with spin-up electrons, while for 
— 1 < P < 0, the polarization is reversed. 

For a specific demonstration, we consider a geometry with the length of the constriction 
L x = 86a, the width of the wide region L y = 70\/3a, L = 34a, and the width of the 
constriction W = 40\/3a; we adopt realistic parameters given in Table I for silicene and 
take A v = E c l for the wide regions to be in the SVPM phase, and A v = 0.053t (> Aso) and 
h = 0.05t for the gated constriction to be in the marginal VPM phase. We set Ep = 0.07t, 
resulting in 2N = 10 right-moving modes. The resulting spin polarization as a function of 
the effective chemical potential, fio = Ep — Uq, is shown in Fig. [2b. Clearly, for positive 
/io < 2h ~ O.leV, the current flows entirely within the conduction band of the left valley 
and the polarization reaches almost 100% efficiency with down-spin. Conversely, for negative 
Uq > —2h, the current flows within the valence band of the right valley and nearly 100% 
spin-polarization is now achieved with up-spin. The case discussed here is the most optimal 
one. It is possible to use different E z and h values in the constriction, which may turn silicene 
into other phases, but we find that the spin filter still remains functional (See Supplementary 
Figs. S3-S5). 

Transport properties. In order to extract which transport features are the key ingre- 
dients for the high-efficiency spin polarization in our proposed device, it is instructive to 
consider a simple transport arrangement which bisects a long strip of silicene (with zigzag 
edges) into two halves as shown in Fig. [3h- We now apply out-of-plane electric fields E\ 
and E2 to the left and right regions, respectively. For simplicity, we assume the Zeeman 
field h = hereafter. The robust helical edge states (which cross from the conduction band 
to the valence band in the QSH phase as Ei < E c ) may or may not be present depending 
on the strength of the fields. Let us examine the case, < E c < E x with E 2 = —E\ (See 



Fig. |3b). Making E 2 negative but still equal to E± does not change the energy dispersion, 
but switches the spin configurations around the two valleys. Our calculation of the con- 
ductance and spin polarization (See Supplementary Table SI) indicates that for one specific 
spin polarization (e.g., up-spin), intervalley scattering is still partially allowed, although 
with serious suppression of the total conductance. But, spin-flip scattering becomes entirely 
forbidden. This is in sharp contrast to the classic example of the valley filter in graphene 
where valley index is preserved during the scattering process due to special features of the 



well-known single- valley edge mode 23|, |33|, |34j . Remarkably, in our proposed spin filter the 



bulk transport from the coupled valley and spin degrees of freedom is the key instead of 
the edge states. In other words, bulk states with non-spin flip scattering processes are the 
crucial ingredients for the high efficiency spin polarization in our case. 



III. DISCUSSION 

Our high efficiency, field tunable spin filter based on silicene takes advantage of charge 
carriers in the bulk system, small Rashba SOC, and controllable spin splitting due to IS 
breaking, A v . We expect therefore our spin-filter to be robust against weak disorder (com- 
pared to Aso + h) and edge imperfections. To test this expectation, we have added random 
onsite potential with strength less than 0. It (~ h) to the constriction. The polarization dia- 
gram is found to be almost quantitatively the same as that shown in Fig. |2b. In particular, 
suppression of the polarization is small, less than 2% as | /xo | < 2/i. Furthermore, we have also 
changed the connection area of the constriction to the wide regions to be armchair-like, and 
introduced random vacancies around the edge of the constriction (up to 20% of edge atoms 
removed), and found that the degradation of polarization is again not severe (less than 1%) 
for such perturbations (See Supplementary Figs. S5-S7 and Supplementary Discussions). 

The working temperature of our proposed silicene spin filter will be controlled by the 
parameter 2Aso = 8.4 meV, yielding a temperature of 97K which lies above the boiling 
point of liquid nitrogen. A similar germanene based device with 2Aso — 23.6 meV could 
be operated at room temperature. As to the characteristic gate voltage, we used a fairly 
large value of 1.57 V/A for E z in the present simulations corresponding to h m 0.052 eV in 
the narrow constriction. Such a sizable exchange field could nevertheless be induced by the 



magnetic proximity effect as shown in the case of graphene. 



35] Using a few times smaller 



value of E z with a corresponding smaller h (while maintaining the marginal VPM phase) 
will still keep our spin-filter functional, although the efficiency and the working region for 
(jlo could be moderately suppressed due to decrease of spin-splitting in the subbands. We 
emphasize however that such reduction in efficiency could be minimized by elongating the 
constriction, as indicated by the simulations of Supplementary Table S2. 

Our spin filter device derives its unusual transport properties via the 2D features of 
silicene under IS breaking and interplay between the electric and magnetic fields. To further 
elucidate the uniqueness of silicene and its potential for future spintronics applications, we 
discuss the Y-shaped spin/valley separator shown in Fig. [H The purpose of the device is to 
separate the two spin/valley polarizations from the incoming lead 1, with one running to lead 
2 and the other running to lead 3. This separator could be operated as follows: First turn on 
an out-of-plane electric field, E z > E c , in the central silicene sheet, tune chemical potential 
\i into conduction bands, and then apply an in-plane electric field by setting potentials, e.g., 
V\ > Vi = V3, at the terminals of silicene. This setup causes charge carriers to acquire an 
anomalous velocity proportional to the Berry curvature in the transverse direction, similar 



to that found out by Xiao et al. in graphene 36] . By linear response theory with negligible 



Ar, the Hall conductivity around valley K,, (rj = ±1 for two valleys) with spin polarization 
a is given by 

< = 4h {l - K ~1° XS ° )- (5) 
In other words, such a setup leads to valley and hence spin polarization imbalance at output 
terminals V2 and V3 (with opposite polarization between them), and results in non- vanishing 
valley Hall [<xJj Vallcy) = ^(1 - f)} and spin Hall conductivity [a^ pln) = ^f(^f)], respec- 
tively. Detailed demonstration of this separator requires delicate simulations, which will 
be considered elsewhere. Our study demonstrates that silicene (and related IVA group el- 



ements with honeycomb structure [17] ) provides a great potential host for manipulating 
spin / valley degrees of freedom efficiently, moving us a step closer to realizing the dream of 
spintronic / valleytronic applications. 

IV. METHODS 

Our first-principles calculations are based on the generalized gradient approximation 



(GGA) 



371 ] using full-potential projected augmented wave method 



381 ] as implemented in 



the VASP package 39]. The 2D low-buckled honeycomb structures of Si, Ge, Sn, and Pb 
were optimized using a 30 x 30 x 1 Monkhorst-Pack k-point mesh over the Brillouin zone 
(BZ) with 350 eV cutoff energy. We note that our computed value of band gap of 8.4 meV 
in silicene is close to the value of 7.9 meV reportedjn Ref. 
are based on the iterative Green's function method 
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The transport simulations 



32| (See Supplementary Methods). 
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TABLE I. Model parameters. Fitting results from first-principles calculations for stand-alone, 
2D Si, Ge, Sn, and Pb sheets. Relevant parameters are lattice constant (a) with buckled structure, 
buckling distance (A c ), transfer energy (t), spin-orbit coupling (Aso)> Rashba spin-orbit coupling 
(Ar), and the linear dependence of the applied electric field on A v (I = X v /E z ). 

a A c t A S o A R I 
(A) (A) (eV) (meV) (meV) (e-A) 
Si 3.87 0.44 1.04 4.2 8.66 0.035 
Ge 4.06 0.68 0.97 11.8 2.81 0.046 
Sn 4.67 0.84 0.76 36.0 18.75 0.055 
Pb 4.93 0.89 0.72 207.3 0.06 0.143 



FIG. 1. Nearly fully spin-polarized states of silicene. a, Low-buckled 2D honeycomb struc- 
ture of silicene. Due to the buckling, the inversion symmetry can be removed by an external 
out-of-plane electric field E z or when the thin film is placed on a substrate, b, Schematic spin- 
resolved band structure of silicene around K and K' points in the presence of an out-of-plane electric 
field, E z < E c . The red and blue arrows indicate the spin direction. The band structure (c) and 
spin-resolved wavefunction (d) in the critical phase (E z = E c ). The tight-binding model (green 
lines) is seen to faithfully reproduce the band dispersion and the nearly 100% electron intrinsic 
spin-polarization of states near the K-point obtained by first-principles calculations (red dots), e, 
The spin splitting energy for silicene, and Ge, Sn, Pb counterparts as a function of E z obtained by 
first-principles calculations. 

FIG. 2. Tunable high efficiency spin filter, a, Geometry of the spin filter and the profile of 
the potential barrier U{xi). The two colored atoms on the lattice emphasize the buckled structure, 
b, Spin polarization of the filter as a function of fiQ in the constriction. The blue (green) line 
corresponds to the case of applying potential barrier with a rectangular (smooth) shape (see text), 
c and d are typical dispersion relations for the wide and the constriction regions, respectively. 
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FIG. 3. Transport selectivity due to spin and valley degrees of freedom, a, Transport be- 
tween two separated regions (by applying different electric fields) of a silicene strip, b, A schematic 
illustration for robust spin-preserved (instead of valley-preserved) scattering processes. 



FIG. 4. Y-shaped silicene spin separator. A schematic Y-shape separator of silicene with 
three local gates (leads). The current separates into the gate 2 and gate 3, respectively, carrying 
opposite spin/valley degrees of freedom. 
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